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Context. Low-mass satellites, like asteroids and comets, are expected to be present around the black hole at the Galactic 
center. We consider small bodies orbiting a black hole, and we study the evolution of their orbits due to tidal interaction 
with the black hole. 

Aims. In this paper we investigate the consequences of the existence of plunging orbits when a black hole is present. 
We are interested in finding the conditions that exist when capture occurs. 

Methods. Earlier analysis of the evolution of classical Keplerian orbits was extended to relativistic orbits around a 
Schwarzschild black hole. 

Results. The main difference between the Keplerian and black hole cases is in the existence of plunging orbits. Orbital 
evolution, leading from bound to plunging orbits, goes through a "final" unstable circular orbit. On this orbit, tidal 
energy is released on a characteristic black hole timescale. 

Conclusions. This process may be relevant for explaining how small, compact clumps of material can be brought onto 
plunging orbits, where they may produce individual short duration accretion events. The available energy and the 
characteristic timescale are consistent with energy released and the timescale typical of Galactic flares. 

Key words. Galaxy: nucleus - Galaxies: active - Physical data and processes: black hole physics 



1. Introduction 

Most galaxies, if not all, harbor massive black holes at their 
center. Some of them manifest their presence in a violent 
way and show up as active galactic nuclei; others, like the 
one in our Galaxy, are quiescent. The mass of the black 
hole at the center of the Milky Way has been estimated 



as Mbh = j.di =e u.jz x iwivifT) irom oDservauons wren 
SINFONI ( Eisenhauer et al.|2005 1 confined within a radius 
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3.61 ± 0.32 x 10 b Mo from observations with 



of 45 AU. The Galactic center is only w 8 kpc from the 
Sun and therefore it is the nearest laboratory where we 
can study the environment of super-massive black holes in 
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(2007 



The standard assumption for AGNs is that gas accretes 
onto the black hole forming an accretion disc and releasing 
energy. On the other hand, some non-active ga laxies exhibit 
short-timescale, low-intensity flaring activity ( Genzel et al. 
2003 |Belanger et al.| [20061 |Eckart et aL"1|2006| |Meyer et 



al.||2006 ), which could be the result of individual accretion 



events of small bodies onto the central black hole. 

It seems reasonable to assume that stars at the Galactic 
center are surrounded by planets and by other small orbit- 
ing bodies, like asteroids and comets. Such low-mass satel- 
lites (LMS) may be stripped off their parent stars by tidal 
interaction, while they approach the black hole. In this way 
they contribute to a distribution of low-mass objects that 
cluster the central black hole. According to |Cochran et al.| 
( 1995 1, the Edgeworth-Kuiper belt of our Solar System may 



Send offprint requests 
andrej . cadezOfmf .uni-lj .si 



to: 



A. Cadez, e-mail: 



still contain as many as 2 x 10 8 objects with radii < 10 km. 
One might therefore expect that LMS will be copiously 
present all the way down to the black hole. 

In this paper we study the evolution of orbits of LMS 
caused by tidal interaction with the black hole. We extend 



the analysis presented by Hut (1980 1981 1982) concern 



ing the tidal evolution of close binary systems in the weak 
friction model in the Newtonian approximation. We show 
that orbital evolution of these objects leads to capturing or- 
bits by the black hole and argue that tidal interaction can 
rapidly inject enough energy to explain the energy source 
of Galactic flares. 



2. Low-mass satellites in the Galactic center 

2.1. Energy-loss mechanisms 

The LMS are injected into the asteroid cloud with low rel- 
ative velocities with respect to the donor star, therefore 
their average random velocity is roughly the same as that 
of stars. In thermodynamic equilibrium, LMS would be ex- 
pected to reach much higher random velocities; however, 
it is not difficult to show that such equilibrium cannot be 
reached, since the relaxation time i a * to exchange energy 
between stars and LMS is 
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where er a is the random velocity of LMS, n*, <7*, and r g » 
are the density, random velocity, and gravitational radius 
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of stars in the central cluster, A is the ratio of the maxi- 
mum and minimum impact parameters, and Xh is the non- 



resonant relaxation time from Hopman & Alexander ( 2006 



hereafter HA06). Since Th is already ~ 10 M yr, a a cannot 
increase much before the relaxation time is longer than the 
Hubble time. This means that LMS do not reach thermody- 
namic equilibrium random velocities, and the Bondi radius 
of the LMS mass distribution is approximately the same 
as that of the stellar distribution, i.e. re ~ 2 pc (HA06). 
However, the direction of angular momentum of such LMS 
changes on a timescale shorter than Th, therefore the dis- 
tribution of these objects eventually becomes spherically 
symmetric around the central black hole. 

Just like stars, the LMS eventually accrete onto the cen- 
tral black hole. To do so, they must lose orbital energy 
and/or angular momentum. While stars eventually lose or- 
bital energy by emitting gravitational waves, LMS are not 
massive enough for this mechanism to play an important 
role. Cons idering the gravitational radiatio n energy loss 
timescale ( |Misner, Thome fc Wheeler ||1971[ Eq. 36.17b) 
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> 10 1 



yr 



10^g 



(2) 



where M and r g (= GMbh/c 2 ) are the mass and gravita- 
tional radius of the black hole, m the mass of the LMS, 
and a(> 4r g ) the radius of its orbit. One finds that satel- 
lites with mass less then 10 23 g can circle the central black 
hole even at a — 4r g for the whole Hubble time. 

However, hydrodynamic drag by circum-black-hole 
cloud might become important in dissipating the energy 
of LMS. T he timescale can be estimated using similar ar- 
guments to Narayan] ( |2000[ ) to obtain 
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where i? e n is the effective radius of LMS, and /3 gas the den- 
sity of the interstellar gas. In the second line we assume 



that the satellite density is 5 g/cm 3 



4 r g , and p g 



expressed in terms of hydrogen atom density. The energy 
loss due to drag strongly depends on the density of the gas 
and becomes comparable with nonresonant relaxation, if 
the mass of LMS is below 



rn = 2.24 10 15 g 



"H 



10 3 cm 3 



4 r„ 



(4) 



For example, an object circling at the Bondi ra dius in a 
gas with density nn ~ 20 cm~ 3 ( |Xu et ah 2006) must be 
less massive than ~ 5 g for hydrodynamic drag to take 
over. In the very vicinity of the black hole, the gas density 
may be higher and also the factor a/r g is much larger, yet, 
since there is no evidence of an accretion disk (Alexander 
|2005[ ) , it does not seem plausible for these factors to increase 
enough for hydrodynamic drag to play a very important 
role. 

Since the evolution of orbital parameters is dominated 
by the process with the shortest timescale, we estimate that 
nonresonant relaxation occurring on a timescale of ~ 10 9 
years (HA06) is one of the most important mechanisms 
until the LMS approaches the black hole to within a few 
Roche radii, where tidal interaction takes over. 



2.2. The population of low-mass satellites 



It has been determined by Paumard et al. (2006) that the 



density of stellar distribution increases as l/r z to within 
the inner 1" of the Galactic center. At this distance, corre- 
sponding to 2 x 10 5 r g , there is a sharp drop, interpreted by 
HA06 as the radius where gravitational radiation extracts 
orbital energy at a sufficient rate to clean the inner region. 
Since for low masses gravitational radiation energy loss is 
negligible, we expect that their density keeps inc reasing ei 



ther as n(r ) cx r (HA06) or as n(r) oc r (Paumard 



et al.||2006[ ) deep down to the very vicinity of the black 
hole, where tidal effects may extract orbital energy and an- 
gular momentum. Note that there are no stars (very few) 
and no gas (no steady emission lines) in the central region 
(r < 2x 10 5 r g ), thus it is likely populated predominantly by 
low-mass satellites. Hence, we assume that the dynamics of 
this central cloud is dominated by the exchange of angular 
momentum by nonresonant relaxation (HA06) and by tidal 
interaction between the object and the central black hole. 
As the nonresonant relaxation timescale is roughly position- 
independent and the tidal interaction grows inversely pro- 
portional to a high power of the distance from the black 
hole, the tidal interaction finally prevails in determining 
the way in which accretion occurs. 

One should expect a fair proportion of LMS to move on 
highly eccentric, low-periastron orbits. Tidal forces do sig- 
nificant work on such satellites near periastron. This lowers 
their orbital energy and starts the significan t evolution of 



orbital parameters (Gomboc & Cadez 
text two classes of such satellites shou 



2005]). In this con- 
d be taken into ac- 



count: those that are gravity dominated (i.e. those whose 
fundamental quadrupole frequency is v g ~ 1\J Gp/3ir) and 
those that are solid-state dominated (whose fundamental 
quadrupole frequency is v s ~ jC s /R). (Here p is the density 
of the body, R its radius and c s the speed of sound.) Taking 
c s w 5 km/s and p ss 5 g cm -3 as typical values, we find 
that the radius dividing the two classes is close to the radius 
of the asteroid Ceres. Therefore all gravity-dominated satel- 
lites should have about the same fundamental quadrupole 
frequency, corresponding to the period of about 54 min- 
utes. All smaller satellites should have shorter fundamen- 
tal periods. This means that gravity-dominated satellites 
start rapid tidal evolution when their periastron reaches 
r p = (GAfbh/^TWg) 2 ) 1 / 3 10r g . Solid-state dominated 
bodies may start significant tidal evolution even closer to 
the black hole. 



3. Tidal evolution of the orbits 

Significant tidal orbital evolution for gravity-dominated 
solid bodies will start when their periastra reach down to 
«9 r g . We show, however, that solid-state dominated bod- 
ies are also strongly affected by tides and are expected to 
be heated by them, so that at a certain stage, they are 
expected to melt and also to become gravity-dominated. 

The first stages of tidal evolution of the orbit can be 
investigated using Hut's formal ism for the spin-orbit evo- 



1980 



1981 19821, at 



lution of the two-body system ([Hut Ki 
least until relativistic regime is reached. The evolution is 
governed by the parameter a, which is the ratio of the 
orbital and rotation angular momentum that the binary 
would have at stable equilibrium, characterized by ao, the 
radius of the stable circular orbit and ujq, the orbital and 
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spin frequency of both bodies. The parameter a is deter- 
mined by the (conserved) value of the total angu l ar mo - 
mentum of the system by Eq. (57) given by |Hut| ( |1981 1. 
We take the s pin of both bodies into account, since it has 
been shown by |Fang fc Love lace (2005J that "the black hole 



absorbs angular momentum and energy at the same rate as 
the moon's tidal field sends energy and angular momentum 
into the hole's horizon" . Thus we generalize Hut's equation 
to 

L = [*(GM)H(m l rf 1)g Ri 



m 2^ 2)g i?|)//x)i(a^ 



a 4 J , 

where mi, R\, f(i) g , and 7712, i?2, r (2)g are the mass, radius, 
and gyration radius of the two bodies, M. — mi + and 
[i = m\m 2 j M.. 

The (conserved) angular momentum, calculated ini- 
tially, is 



1 R 2 cu 2 + n(GMr Q )* , (6) 



■ (1)gJ . H wi -r m 2 r (2)g 

where uj\ and u>2 are the two initial spin angular velocities, 
h = /i(GAiro) 1 / 2 is the initial orbital angular momentum 
with r = a(l — e 2 ), and a and e are the semi-major axis and 
eccentricity of the initial orbit, respectively. Let m\ 3> rri2 
and Ri rm g 3> R 2 ?"(2)g and define 

Z = (m 2 rf 2)g Rt)/(m ir f 1)g Ri) (7) 

Y = h/(m 1 r 2 1)g R 2 u Jl ) . (8) 

Equating Eqs. (|5| and ([6]), we obtain 

(M/mi)5(GM/iJ?)»(l + Z)»(ar*+cr*) 

= lu 1 (1 + Zu; 2 /lu 1 + Y) . (9) 

Note that fi<j = (GM/i? 3 ) 1 / 2 is the angular frequency at 
which m 2 would orbit mi at its gyration radius. If one 
assumes that mi = -Mbh and m 2 = m, then = 
ciKerr/^gC where OKcrr is the Kerr angular momentum pa- 
rameter. Since ^> m, we may replace M = Mbh, 
fi = m and write R 2 — R. Finally, from the above equations 
we obtain 



(a 1 







(1 + Zlu 2 /lu 1 +Y)(1 + Z)-i 



From this expression it is clear that the angular momentum 
of the black hole determines a; moreover, since aKerr/^gC 
is not expected to be exceedingly small and Z and Y are 
very small, the equation for a is essentially 



(a 1 







a Kcrr 

r„c 



\ m 



(11) 



As the righthand side is very large, it follows that 
a is either a very small or a very large num- 



ber. Since, by definition, 



m 2r 2 2)g R 2 )uJo 



= h [(m ir 2 1)g R 2 



|(m/Mbh)(ao/?"g) 2 it is clear that 
we are not interested in very high, but in very low values of 
Hut was not interested in such orbits. We note that, if 

1 The equilibrium orbit is circular with radius ao, so that ho = 
ma^u>o and, since it is co-rotating with the black hole, ujo = 
UKerr- With respect to the moment of inertia of the black hole 
m i r (i) 9 ^i = Jbh/^Kcrr, the moment of inertia of the asteroid 
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Fig. 1. Orbital evolution tracks in the eccentricity (e)- 
periastron distance(r p ) plane. Different types of evolution 
are displayed. Arrows indicate the direction of time. The 
left panel shows tidal evolution toward stable corotating 
circular orbits at r' — 450r g . The right panel is an en- 
largement of the left bottom part. It shows the evolution 
of very low orbital angular momentum orbits. (See text 
for details.) In this example, intended to show the topol- 
ogy of small a tidal evolution, the value of a is arbitrar- 
ily taken to be 0.05. Different curves correspond to ini- 
tial f p = 25, 50, 100, 200, 250, 400, 600, 800, 1000 (left panel) 
and f p = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.85,0.9 (right 
panel) . 



r p (e\a') is the solution to the equation (27) in Hut (19821 
for a' > 3 and if "fr p (e\a) is the solution of the same equa- 
tion for a < 3 such that a' and a are the two solutions of 



Eq. (57) in Hut (19811, then a scaling factor ^exists, such 
that both are solutions of the same equatiorjj Therefore 
both a's solving Eq. (Ill produce the same solutions. 



The equation for the flow of trajectories in the (r p ,e) 
plane for very low values of a is the same as that for very 
high values of a and can be written using Eq. (27) in Hut 
(10) fli982l ) 



dr p _ f p (-49 + 81e - 70e 2 + 22e 3 - 17e 4 + s 5 ) 

~de ~ ll(-2 + e)(-l + e)(21 - 28e + 18e 2 - 4e 3 + e 4 ) ' 

Here r p = r p /ao is the normalized minimum separation of 
the two objects at periastron, f p = a(l — e)/ao- Noting that 
e = 1 — e, where e is the eccentricity of the orbit, the above 
may be integrated to 

f p =(21 - 28(1 - e) + 18(1 - e) 2 - 4(1 - e) 3 + (1 - e) 4 )n 



x e 1 * 1 (1 + e) 1 f n 



(13) 



where r n is the integration constant. The meaning of these 
limiting solutions for very small a is i llustrated i n Fig. [T] 



by the numerical solution of Eq. (27) in Hut ( 1982~| . In this 
example, where the value of a was chosen equal to 0.05 



can be neglected. According to Ashtekar fc Krishnan | ( 2004 1 , 2 The exchange a — > a' changes the unstable equilibrium 
Jbh/wKcrr = -A^bh^h! where i?h is the horizon radius of trie black state f_ to the attractor f p = 1 and vice versa with respect to 
hole, which becomes 2rv in the limit of small axcrr- 



the tracks shown in Hut ( 1982 I 
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Fig. 2. Angular momentum (A = 2r ^ r ) - orbital energy 

(E) during the tidal evolution of the orbit for the evolution 
tracks shown in the right panel of Fig. [I] 

(a' = 9720)]^] the equilibrium radius a' is at f p w 450 
and only orbits with f p < 0.8 shrink as they circularize. 
In Fig. [T] the left panel shows tidal evolution toward sta- 
ble corotating circular orbits at r' — 450r g . These orbits 
are initially quite eccentric and have low spin (the upper 
five curves) . The remaining orbits start as circular but with 
large spin. If the initial spin energy is high enough, it can 
be transferred to orbit, first elongating it and then circu- 
larizing again after joining the circularization track (two 
upper curves starting to the left). If the initial spin en- 
ergy is not high enough, the circularization track cannot 
be reached and the orbit keeps elongating at a slower and 
slower pace (lower two curves). The right panel is an en- 
largement of the left bottom part. It shows the evolution of 
very low orbital angular momentum orbits. Very high ini- 
tial spin angular momentum and energy can be transferred 
to orbital momentum and energy, leading to higher perias- 
tron orbits (the upper three curves). In the case of small 
initial spin, the orbital evolution starts with orbital energy 
dissipation by tidal interaction and little angular momen- 
tum transfer, leading to less eccentric orbits with shorter 
and shorter orbital periods. As a result, the orbital angu- 
lar momentum eventually transfers to spin as the object is 
forced into faster and faster corotation. As a' — > oo, the 
corresponding a' also tends to infinity and only the shrink- 
ing orbits remain as candidates to bring the body down to 
the horizon of the black hole. With the solution for r p and 
noting that 

h = m(GM hh r ) * = mr g c((l + e)r p /r g ) 5 (14) 
E = -GM hh m(l - e)/2r p , (15) 

one can calculate tracks of orbital evolution in the angular 
momentum - energy plane (see Fig. |2j. Using these tracks, 
one can calculate the evolution of orbital radial turning 
points. 

We do this using both the Keplcrian and the black hole 
effective potentials. For simplicity we only consider the 
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Fig. 3. The evolution of radial turning points in the Kepler 
(red) and Schwarzschild effective potentials (black) for the 
case of a highly eccentric orbit starting with periastron at 
9.5r g and apoastron at 110r g . The right figure is an enlarge- 
ment of the final evolution. The gray line indicates the third 
root of the effe ctive potential equation (Misner, Thorne & 
Wheeler 



1971) 



Schwarzschild metric (Fig. [3|. The important qualitative 
difference between the two cases occurs for low values of or- 
bital energy and angular momentum. In the Keplcrian case 
the effective potential guarantees two turning points for all 
bound orbits, since it has one minimum and no maxima, 
while in the Schwarzschild case the effective potential has 
a minimum (V m i n — V (r m j n )) and also a maximum (V max ) 
at small r max < r min (Misner, Thorne & Wheeler 1119711). 



This is an unreasonably large number with respect to our 
problem, but sufficiently small and convenient to show the topol- 
ogy of the solution space. 



Therefore, the inner turning point disappears in a paramet 
ric family of orbits, when the effective potential maximum 
becomes less than the orbital energy. For the parametric 
family of orbits shown in Fig. [2] this happens when r max is 
just slightly larger than 4r g . In Fig. [3] we then show that, 
from the point of view of turning points, Schwarzschild or- 
bits do not necessarily circularize, but the inner turning 
point disappears when the outer one may still be at ss 20r g . 
The relativistic theory of tidal orbital evolution may pro- 
duce somewhat different tracks in the A — E plane, yet we 
believe it is reasonable to expect the results not to change 
qualitatively, because the exchange between orbital and in- 
ternal angular momentum and energy occurs very locally 
in the small volume of the tidally distorted body. In any 
case, when the total energy of the body closely approaches 
the maximum value of the effective potential (Vm ax ), the 
body experiences so shallow an effective potential that it 
does not return to the apoastron, but winds about r max . 



4. Time scales for tidal evolution 

To estimate the timescale for tidal evolution of the orbits, 
consider the energy loss per orbit as given by Eq. (A10) of 
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Hut (1981 4 It can be written as 
{GmM hh R 2 ' 



AE = - 



9tt 



-fee 2 



'GM hu (r R 



(16) 



where r R = [M^/m^'^R is the Roche radius. The non- 
resonant energy timescale is then 



which becomes 



^EN — — iorb — /AE 

la 



64(l + e)^r8 



(17) 



9c 2 y/T^e e 2 (64 + 240 e 2 + 120 e 4 + 5 e 6 )kR 2 r g r R T ' 

(18) 

where i or b is the orbital period, k = 0.75 the apsidal motion 
constant of the primary (considered as an incom pressible 



19811. 



fluid) , and r a constant small time lag ( Hut 

Remember, however, that Hut only considers weak tides 
acting on the deformed body m. Such tides are quite off 
resonance, and this would not be the case in the proximity 
of the black hole, where tides would certainly be resonant if 
the body liquefies. Therefore a better es timate is as follows. 
For re sonant tides the body liquefies. Gomboc fc Cadez 
( 2005 1 obtain the following equation for tidal energy per 



penastron passage, valid in the limit of high eccentricity, 
/GM hh mR 2 



AE 



(19) 



with the resonant timescale for tidal evolution, 

..2 



(l-e)r 2 
'2R 2 e 2 (/3) 



(20) 



and e 2 ((3) is a function that peaks at resonance, i.e. at (3 ~ 
1. Here (3 = r R /r p is the Roche penetration parameter. For 
easier comparison we rewrite Eqs. 18 and 20 into a simpler 
form 



t E N = 1-4 10 11 yr x 



Qfm{e) fW km 



ER 



2.9 10 10 



yr x 



e 2 (p)(3i 



R 

10 km 
R 



(21) 
(22) 



where Q is the resonant damping factor, r = (loqQ) , 
loq = \J GMbh/^R, and po — lg cm -3 . The function /en(c) 



is 



/en(c) 



16(1 



25 v / T Tr e e 2 (64 + 240 e 2 + 120 e 4 + 5 e 6 ) ' 
and /er(g) (valid only for highly eccentric orbits) is 

/er ( s) 



(23) 



(24) 



Both functions are of order unity on the interval [0.1,0.9]. 
In Fig. [4] we plot £en and £er as a function of f3 for R = 
10 km, Q — 100 and p — po, assuming that the body is 
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Fig. 4. Resonant and non-resonant timescales as a functio n 
of $ for p = po, R= 10km, and Q = 100 (see Eqs. |21|22 1. 



liquefied, 



gravity-dominated. For example, if an R 
g cm 

<er = 3 x 10 10 yr x (10/30) 2 (l/5) 7 / 6 = 5.1 x 
this case <en is approximately twice as long. 



30 km object with p = 5 g cm 3 is on a f3 = 2 orbit, then 



10 s 



yr; 



f(jr 



5. The fate of a small body 

Consider now a cold solid object that is scattered on an 
orbit with a periastron of a few 10 r g . Tidal energy loss 
and the corresponding energy timescale can be accurately 
calculated using Hut's formalism, since tides on such an 
object would be well below resonance. However, one must 
replace the apsidal motion constant k with a lower value 
and loq with the higher angular frequency of the quadrupole 
mode of the solid body. Thus, to obtain t he e nergy time 
for such a solid body, the timescale in Eq. (18 1 should be 
multiplied by the ratio (^ s /i/ g ) , to obtain 



^ES 



= 5.7 10" yr x ( ^ V 



1? 



The average heating power is 



-Pheat — \AE\/t OI b 
1 

2 



(25) 



(26) 



-GMbh tides' 2 



The typical 
T t' _1 , where D 



where the last equality follows from Eq. (|17l 
thermal diffusion timescale is T so iid = R 2 ir 
is the thermal diffusion coefficient, which is on the order 
of 10~ 6 m 2 s _1 for most rocks, so that r so iid - 3 x 10 5 
years for a 10 km object. It is a simple exercise to calculate 
the central temperature that such an object would reach in 
equilibrium when the heating power is thermally diffused 
to the surface, which is black-body radiating. Neglecting 
the surface temperature, one obtains T c = Pi lcat (87ri?A c ) _1 , 



4 For the sake of simplicity we neglect the contribution of work 
done by the tangential component of tidal force, which depends 
on the difference between the orbital and spin period. 



(i 
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where the heat conductivity A c ss 1 — 2Wm _1 K _1 is typical 
of most rocks. Thus rocks could reach central temperatures 
on the order of 



T c = Q.Q2K 



1 



Q/en(6) 



10 k m ; \poj 



(27) 



Assuming that natural values for Q are a few hundred up to 
a few thousand, the same as for normal modes of Earth, and 
noting the extremely steep dependence of T c on R, (3, and p, 
one expects a few 10 kilometer rocks to melt (T c > 3000-ftT) 
at their centers a few hundred thousand years after entering 
a (3 > 2 orbit. Melting weakens the body, increases its cou- 
pling to tides, and intensifies heating. Therefore, the core 
melts farther out, but does not essentially increase its tem- 
perature, because the convection that sets in is extremely 
efficient in transporting heat from the center to the remain- 
ing solid envelope. This runaway process stops when further 
melting can no longer weaken the body, i.e. when the body 
is no longer solid state but gravity-dominated. If melting 
had started at /? > 1, then the body finds itself at tidal 
resonance, and the pertinen t ti mescale shortens consider- 
ably and becomes Ter (Eq. 22 1. The surface temperature 
of the body is now determineaby the balance between tidal 
heating and black body cooling, so it can be expressed as 



bb 



610 K 



(l-e)e 2 (/?) 
/er(c) 



R 



10 km 



(28) 

Because of the 1/4 power in Stefan's law, the black body 
surface temperature is a slowly varying function of R, (3, 
and p. Roughly speaking, it is expected to be just below 
the boiling temperature of rocky material (a few thousand 
Kelvin) for melted rocks at a tidal resonance whose size is 
between about 10 and 100 kilometers. Much smaller rocks 
would not melt at all, and much larger rocks, moons, and 
planets would melt farther away from the black hole and 
would evaporate or split into smaller pieces or drops when 
reaching the Roche radius. Evidence of a si milar tidal melt- 
ing m ay also be found in our Solar System (Dermott et al 
19881) . 

Heating power increases slowly with slowly increasing 
f3, and the melted body undulates as it moves between pe- 
riastron and apoastron on an eccentric enough orbit that 
keeps it outside the Roche radius most of the time. 

The above considerations apply until angular momen- 
tum transfer becomes important. Angular momentum can 
be transferred from orbital to spin of the black hole and 
from orbital to spin of the body. The angular momentum 
transfer to a Schwarzsch ild black hole has been studied by 
Fang & Lovelace (2005) for the case of tidal coupling of a 
black hole and a circularly orbiting moon. They find that 
tidal interaction gives rise to orbital angular momentum 
loss rate (see their Eq. 27): 



32 m 



5 M bh r« 



•-he 



(29) 



In our case m/M^ is so small that the timescale result- 
ing from Eq. (29 1 is longer than the Hubble time, so this 
effect is negligible. On the other hand, the angular momen- 
tum transfer to the spin of the body can work. We are 
now considering an already melted electrically conducting 
object that is constantly elongated by the mounting tidal 



force. Therefore, its moment of inertia is increasing, while 
the ever longer object is spinning with the orbital angular 
velocity. This process only has meaning if the body remains 
a whole and sufficiently rigid, i.e. if there is a mechanism 
that distributes angular momentum to all parts of the body. 
In our case sufficient rigidity can be provided by a frozen-in 
magnetic field. The necessary field must be strong enough 
to insure that Alfven waves can travel the length (I) of 
the body before the frozen-in magnetic field decays, i.e. 
«ATm > I- Here va ~ B/y/pop is the Alfven velocity and 
T m = ct/j,o^ is the magnetic field decay time, where l t is the 
transverse dimension of the body. Taking a typical value 



a ~ 4 x 10 6 Q- 



for the conductivity, one finds that 



a field of ~ 100 pG is enough to give rigidity to a 100 km 
object. S uch and much stronge r fields have been found in 
asteroids ( Kivelson et a l.||TT) !),'!). so that it is not unreason- 
able to assume that the LMS in the vicinity of the Galactic 
black hole also possess it. 

Slow orbital angular momentum decay describing the 
lower part of tracks in Fig (2] has no dramatic consequences 
until fully relativistic regime is reached, i.e. until the effec- 
tive potential maximum starts rapidly dropping with loss 
of angular momentum. This happens when the inner turn- 
ing point approaches r max rs 4r g . The shallow potential 
no longer returns the body from periastron to apoastron, 
but the orbit just keeps winding about r max while the con- 
stant large tidal tensor, now experienced by the body, keeps 
squeezing and pulling it apart into a long thread along the 
orbit, exponentially increasing its internal energy. 

We note that the above analysis is an approximation 
somewhat limited by neglecting dynamical effects and rel- 
ativistic aspects of tides. As pointed out in Sect. [3] in the 
relativistic (Schwarzschild) regime, the orbit has no inner 
turning point when the effective potential maximum be- 
comes less than the orbital energy, leading therefore to cap- 
ture orbits with orbital energy that may be considerably 
higher than the orbital energy of the marginally stable or- 
bit. This energy difference is available to do tidal work on 
the body before entering the final turn down the horizon of 
the black hole. Since tidal interaction, which transfers en- 
ergy and angular momentum between spin and orbit, is tak- 
ing place in a very small region of space-time, it might not 
be unreasonable to expect that a full relativistic treatment 
of the problem will not qualitatively change this result, al- 
though it is clear that a number of detailed questions still 
need closer attention. In particular, Hut's analysis should 
be extended to a relativisti c regime above resonance, and 
the |Gomboc fc Cadez ( |2005[ ) analysis would have to include 
higher order modes, mode splitting due to rotation, and a 
relativistic regime. 

6. Conclusions 

In this paper we have investigated the fate of small bod- 
ies, like asteroids and comets, that may find themselves on 
highly eccentric orbits in the vicinity of the Galactic black 
hole. Extending Hut's analysis of the tidal evolution of close 
binary systems, we find that they experience very high tides 
at their periastra, which heat and eventually liquefy them. 
Those objects that are electrically conductive are likely to 
have enough magnetic rigidity to efficiently transfer angu- 
lar momentum from orbit to spin, so that the loss of orbital 
angular momentum leads to smooth transition from bound 
to plunging orbits. We propose this as a mechanism that 



Cadez et al.: On tidal evolution of orbits around black holes 



7 



could bring relatively small and condensed clumps of ma- 
terial to the vicinity of the black hole. Such clumps can 
produce individual accretion events characterized by the 
time scale of the last circular orbit. The tidal energy re- 
leased during this process can reach up to 0.1 mc 2 . This 
means that both the available energy and the characteris- 
tic timescale are consistent with the energy and timescale 
characteristic of Galactic flares. 
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